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ABSTRACT 

We present Spitzer IRAC and IRS observations of the dark cloud LI 287. 
The mid-infrared (MIR) IRAC images show deeply embedded infrared sources 
in the vicinity of the FU Orionis objects RNO IB and RNO 1C suggesting their 
association with a small young stellar cluster. For the first time we resolve the 
MIR point source associated with IRAS 00338+6312 which is a deeply embed- 
ded intermediate-mass protostar driving a known molecular outflow. The IRAC 
colors of all objects are consistent with young stars ranging from deeply embed- 
ded Class 0/1 sources to Class II objects, part of which appear to be locally 
reddened. The two IRS spectra show strong absorption bands by ices and dust 
particles, confirming that the circumstellar environment around RNO 1B/1C has 
a high optical depth. Additional hydrogen emission lines from pure rotational 
transitions are superimposed on the spectra. Given the outflow direction, we 
attribute these emission lines to shocked gas in the molecular outflow powered 
by IRAS 00338+6312. The derived shock temperatures are in agreement with 
high velocity C-type shocks. 

Subject headings: Stars: Pre-main sequence - Stars: Formation - Stars: circum- 
stellar matter - Infrared: stars - ISM: jet and outflows - Individual stars: RNO 
IB, RNO 1C, IRAS 00338+6312 
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1. Introduction 

The dark cloud L1287 contains the galactic nebulae GN 00.34.0 (associated with the 
young F-type star RNO 1) and GN 00.33 .9 which harbors the IRAS point source IRAS 
00338+6312 (d = 800 pc, IPersi et al.l [l988h . Two young stars, RNO IB and RNO 1C0, lie 
slightly south-west of the catalog po sition of the IRAS source . Those objects show properties 
of FU Orionis objects (FUORs): (1) Istaude fc Ne"ckel fjl99lh found that RNO IB brightened 
by at least 3 mag over a period of 12 years and that it showed a variable and blueshifted P 
Cygni profile i n Ha and add i tional broad and double-peaked absorption lines in its optical 
spectrum; (2) iKenyon et al.l (119931 ) obtained near-infrared spectroscopic data for both ob- 
jects and found the FUOR-typical strong 2.3 /im CO absorption bands. As the error ellipse 
of the IRAS source includes both FUORs it was long debated whether there is still another 
deeply embedded source close to the FUOR objects or whether IRAS just measured the inte- 
grate d flux from these two objects . From high-resolution near-infrared (NIR) polarimetric 
maps I Weintraub fc Kastnerl (119931 ) concluded that an additional embedded object should be 



present close to the locati on of the IRAS sourc e. This idea was sup ported by th e discovery of 
a 3.6 cm continuum peak (lAnglada et al.lll994j ) and an H2O maser (lFiebiglll995l ). both nearly 



coinci ding with the IRAS posi tion. A bipo l ar ou tflow in the region was found by ISnell et al. 



(119901 ) and later confirmed by I Yang et al.l (11991I ). As the positions of the IRAS source and 
the FUORs line up along the outflow axis it has long been uncertain which source was 
driving the outflow. From interferometric observations in CS Yang et al.jll995[) s uggested 



that the IRAS source was most likely the driving source. Recently, 
the outflow in CO and came to the same conclusion. However 



Xu et al. 



(2006) m apped 



McMuldroch et al.l (119951 ) 



presented (sub-)millimeter observations favoring RNO 1C as the outflow driving source. 

In this paper we present mid-infrared (MIR) imaging and spectroscopy data taken with 
the IRAC and IRS instruments onboard the Spitzer Space Telescope. For the first time we 
resolve the MIR point source associated with the IRAS source and find additional, partly 
deeply embedded, objects. The two IRS spectra probe the composition of the dense circum- 
stellar environment in the vicinity of RNO 1B/1C and bear additional traces of shocked H 2 
gas. 



We use the same nomenclature as Istaude fc Neckell fjl99lh 
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2. Observations and Data Reduction 

All data were part of the GTO program by R. Gehrz and publicly available from the 
Spitzer data archive. An overview of the observational setup and the datasets is provided 
in Table [TJ The IRAC images were obtained in sub-array mode leading to an effective 
field-of-view of 40" centered on the position given in Table [TJ The spectra taken with 
IRS cover the wavelength range 5 — 37 /im. Overplotting the spectral slits of IRS on the 
2MASS Ks-filter image reveals that apparently RNO IB and RNO 1C were not centered 
within the slits (Fig. [TJ). The short low-resolution spectrum (5.2 — 14.5 /im, R ~ 64 — 128) 
and the short high-resolution spectrum (9.9 — 19.6 /im, R ~ 600) close to RNO IB seems 
to probe mainly flux coming from in between the two objects. The spectrum close to RNO 
1C contains presumably also less flux than expected due to the slight mispointing. The long 
high-resolution part in either spectrum (18.7 — 37.2 /im, R ~ 600) includes flux from both 
components and also additional flux from the IRAS source as all three objects lie within the 
slit of the spectrograph. 

The IRAC images were reduced with the MOPEX package provided by the Spitzer Sci- 
ence Center. Interpolation, outlier detection and co-addition of the images were carried out 
for each filter individually. The astrometry was refined by comparing the position of detected 
sources to known 2MASS objects. The photometry was carried out with the daophot pack- 
age provided within the IRAFg environment. As there was only a limited number of sources 
available within the small field covered by the camera we measured the point-spread function 
(PSF) of RNO 1C and used it as reference to do PSF photometry for all the other sources. 
Following the IRAC Data Handbook 5.1.1 we converted the pixel values from MJysr -1 to 
DNs -1 and computed the corresponding magnitudes via m = — 2.5 log(x) + Azp with x 
denoting the flux measured in DNs -1 and Azp being the zero point for each filteiQ For the 
initial PSF fit we used a PSF size of 2 pixel and applied aperture corrections as described 
in the Data Handbook to obtain the final magnitudes. 

Our final IRS spectra are based on the intermediate droopres (for the low-resolution 
data) and rsc (for the high-resolution data) products processed through the S 13.2.0 version 
of the Spitzer data pipeline. Partially derived from the SMART software package (jHigdon et al. 



20041 ) these intermediate data products were further processed using spectral extraction tools 



developed by the "Formation and Evolution of Planetary Systems" (FEPS) Spitzer science 



2 http:/ /iraf.noao.edu/ 

3 Zero points taken from lHartmann et al.l |2005l) : 19.66 (3.6 fim), 18.94 (4.5 fim), 16.88 (5.8 /im), 17.39 



(8 fim). 
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legacy team (see, Explanatory Supplement v. 3.0 within the FEPS data deliveries^. 

For the short low-resolution observations, the spectra where extracted using a 6.0 pixel 
fixed-width aperture in the spatial dimension resulting in a 39.96 arcsec 2 extraction aperture 
on the sky. The background was subtracted using associated pairs of imaged spectra from 
the two nodded positions along the slit. This process also subtracts stray light contamination 
from the peak- up apertures and adjusts pixels with anomalous dark current relative to the 
reference dark frames. Pixels flagged by the Spitzer data pipeline as being "bad" were 
replaced with a value interpolated from an 8 pixel perimeter surrounding the errant pixel. 
The high-resolution spectra were extracted with the full aperture size (53.11 arcsec 2 for 
the short-high and 247.53 arcsec 2 for the long-high module). The sky contribution was 
estimated by fitting a continuum to the data. Both, the low- and high-resolution spectra, 
are calibrated using a spectral response function derived from IRS spectra appropriate for 
the two nod positions in the slit at which we extract the specta, and Cohen or MARCS 
stellar models for a suite of calibrators provided by the Spitzer Science Centre. The multiple 
orders in our spectra match to within 10%. These small flux offsets are likely related to 
small pointing offsets. We computed correction factors for possible flux loss due to telescope 
mispointing based on the PSF of the IRS instrument. However, accurate results can only be 
obtained if the flux is dominated by one source. The offsets also result in low level fringing at 
wavlengths longer than 20 micron in the low-resolution spectra and at all wavelengths in the 
high-resolution spectra. We removed these fringes using the irsf ringe package developed 
by F. Lahuis. The relative errors between spectral points within one order are dominated by 
the noise on each individual point and not by the calibration. Based on our experience with 
Spitzer IRS data we estimate a relative flux calibration error across a spectral order of ~ 5 % 
and an absolute calibration error between orders/modules of ~ 10 %. These values, however, 
are based on point source observations with accurate telescope pointing. The data presented 
here are more complicated as at least part of the emission is moderately extended, multiple 
sources contribute to the observed fluxes, each contribution is wavelength dependent, and 
the telescope slits were not centered on the main objects. Thus, an accurate flux calibration 
between the different modules is difficult and apparent offsets are further discussed in section 
3.2. 



4 http:/ /ssc. spitzer. caltech.edu/legacy/fepshistory.html 
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3. Results 

3.1. IRAC photometry 

Within the limited field-of-view of the IRAC subarray mode we identified 8 sources 
that were detected in at least three of the four IRAC bands. Figure [2] provides a comparison 
between the RNO 1B/1C region as seen in the 2MASS Ks-filter and the 5.8 /im IRAC filter. 
For the first time the MIR point source related to IRAS 00338+6312 is detected in the IRAC 
band. Additional fainter objects are also present, some of which are detected for the first 
time. Fig. [TT1 shows a color composite image of the RNO 1B/1C complex based on three 
IRAC filters. Table [2] lists all objects that were identified in at least three IRAC bands and 
summarizes the derived fluxes. The errors are based on the results from the PSF-photometry. 
Since we used the PSF of RNO 1C as reference PSF the corresponding errors are relatively 
small compared to the other objects. 



IRAS 00338+6312 and RNO 1GH were not detected in the shortest IRAC band at 
3.6 /mi. The 2MASS image in Fig. [2] shows that at least the IRAS source is deeply embedded 
in a dense dusty environment explaining the non-detection at this wavelength. Fig. [3] shows 
two color-color plots ba sed on the four IRAC bands for all objects listed in Table[2J Following 



Hartmann et al.l (120051 ) who analyzed a large sample of pre-main sequence stars in the Taurus 
star-forming region we use these plots to classify the different objects. The colors of the 
reddest objects (IRAS 00338+6312 and RNO 1G) are consistent with very young protostars. 
The colors of the newly discovered objects RNOl IRAC1 and RNOl IRAC3 are consistent 
with Class 0/1 systems although it has to be confirmed that these are indeed nearby objects 
and not highly reddened background so urces. While the colors o f the FUOR object RNO 1C 



fit in the region of Class II objects from lHartmann et al.l (120051 ) in both plots (Fig. [3]), RNO 



IB fits only in one color-color diagram into this regime. In the right plot of Fig. [3] RNO IB 
is too red for a Class II source in the [4.5]- [5.8] color, but not red enough in the [3.6]- [4.5] 
color to be a Class 0/1 object. Finally, RNOl IRAC 2 and RNO IF, which both show up in 
the 2MASS Ks-band image, can also be interpreted as Class II objects. However, the colors 
of RNOl IRAC2, but partly also of RNO IF, are possibly altered due to local extinction 
effects. 



5 T he object we call RNO 1G appears to be identical to the embedded YSO from IWeintraub fc Kastner 
(|l993h . 



- 6- 



3.2. IRS spectroscopy 

In Figure H] we show the Spitzer IRS spectra without any extinction correction. As 
pointed out earlier the spectral slits of the short wavelength modules were apparently not 
directly centered on the objects RNO IB and RNO 1C and thus the measured fluxes cannot 
be attributed to these objects with very high accuracy Also, there seems to be an offset in 
the flux in the high-resolution regime compared to the low-resolution part of the spectrum. 
To see to which extent this effect is caused by different aperture sizes we plotted the spectra 
in intensities (i.e., flux density per solid angle) rather then in Jansky per wavelength. It 
shows, however, that even after this correction significant offsets remain: Around 13 /im 
the difference between the low-resolution spectrum and the short high-resolution spectrum 
amounts to a factor of ~ 1.25 and ~ 1.08 for RNO IB and RNO 1C, respectively, with the 
high-resolution part showing higher intensities. At least for the RNOl B spectrum this 
offset is larger than normally expected from the calibration accuracy. We thus believe that 
part of this offset can be attributed to the different orientations of the spectral slits on 
the sky probing different regions of the extended emission. The offset between the short 
and long wavelength range of the high-resolution data corresponds to factors of ~ 1.5 and 
«2.0 at 20 /im for RNO IB and RNO 1C, respectively. Also here the short high-resolution 
spectrum shows higher intensities. This can be explained by the approximately five times 
smaller aperture in the short high-resolution module. Although this aperture is probing a 
significantly smaller regions on the sky these regions do presumably contribute in total to flux 
in the above mentioned wavelength regime. The large aperture of the long high-resolution 
module on the other hand is certainly also probing regions without any significant flux (see 
also, Fig. [[J. 



3.2.1. Ices and silicates 
At 6 /xm both spectra show clear water ice absorption bands. At slightly lon ger wave- 



lengt h (~6.85 /im) additional absorption possibly arising from CH3OH ice is present (Ivan Dishoeck 



20041 ) . The 10 /im silicate band is also seen in absorption in both spectra (stronger close to 
RNO 1C), although the shape of the feature differs from the absorption feature caused by 
typical interstellar medium silicate grains. To analyze the differences in more detail we fitted 
a continuum to the 10 /zm region of the spectra and computed the optical depth (Fig [5]). In 
both cases the peak of the absorption is slightly shifted towards shorter wavelengths indicat- 
ing a non-ISM like dust composition. Additional absorption (RNO 1C) or possibly additional 
dust emission on top of the absorption feature (RNO IB) is seen at longer wavelengths. At 
15.2 /im CO2 ice is creating a prominent absorption band and at ~18 /im additional silicate 
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absorption seems to be present. All these features indicate the existence of a dense dusty 
and icy environment in which the two FUORs are embedded. A further analysis of these 
features also for a larger sample of FUORs will be presented in an upcoming paper (Quanz 
et al. in preparation). 



3.2.2. Pure rotational H 2 emission 



In addition to the ice and silicate features, H2 emission lines from purely rotational 
quadrupole transitions are present in both spectra (Fig. H|). While in the spectrum close to 
RNO IB all transitions from S(l) to S(7) can be identified, the spectrum close to RNO 1C 
shows only the lines from S(l) to S(5). This can be explained with the higher continuum flux 
close to the S(6) and S(7) line in the latter spectrum and with apparently lower excitation 
temperatures (see below). The lowest transition S(0) near 28.22 fim is not detected in either 
spectrum. This is due to the strongly rising continuum at longer wavelengths. Here, the 
spectral slit contained flux from both RNO objects and also from the deeply embedded IRAS 
source so that the continuum emission and the related flux errors completely dominate a 
possible weak emission line. A detailed analysis is provided in the appendix. 

Keeping in mind the existence of a molecular outflow that is powered by IRAS 00338 +6312 



and directly oriented in the direction of RNO 1C and RNO IB (e.g., iXu et al.ll2006l ). the 
detection of H2 emission lines in both spectra hints towards shock induced emission related 
to the outflow. The spectra consequently bear information about the circumstellar material 
close to RNO 1B/1C and about the outflow coming from the IRAS source. Since we observe 
H 2 lines even in the 10 /im silicate absorption bands the outflow appears to lie in front of the 
dusty environment as otherwise the high extinction (Ay = 9.2 mag and Ay pa 12.0 mag for 
RNO IB and RNOIC, respectively, IStaude &: Neckellll99ll ) would have prevented a detec- 
tion. The measurement of the relative strengths of multiple H2 lines allows an analysis of the 
physical conditions of the shocked material. For this we assume local thermal equilibrium 
(LTE) and optically thin line emission which is supported by the low E i nstein coefficients 
of the involved quadrupole transitions (Table [3]). Following iParmar et al.l (119911 ) the column 
density of an upper energy level N U (J) is then given by 



N U (J) 



4tt I(J) 



cm 



<r T A ul AE ui 

where I (J) denotes the observed line intensity in erg s" 1 cm -2 sr" 1 , AE u i is the energy 
difference of the two states involved in the transition, A u i is the Einstein coefficient for the 
transition and r is the optical depth at the observed wavelength. In Fig. [6] we show Gaussian 
fits to the observed emission lines in the spectrum close to RNO IB and Fig. [7] presents 
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corresponding fits to the lines detected close RNO 1C. Before fitting the emission lines 
we subtracted the underlying continuum which was fitted with a second order polynomial. 
The finally derived integrated line fluxes and column densities are shown in Table [3J We 
corrected the line fluxes for extinction effects using the results of lMathisI Jl99oh and as s uming 



Ay = 3.55 mag. This value represents the extinction found by iStaude fc Neckell (119911 ) 
towards the nearby star RNO 1 and should be a better estimate for the line-of-sight extinction 
than the above mentioned extinction values towards RNO IB and RNO 1C. In any case the 
influence of the extinction on the derived shock temperatures and column densities (see 
below) is negligible. The 1-sigma errors in the line flux, the line intensity and the column 
density (Table [3]) were derived from generating 500 spectra where we added a Gaussian 
noise to each measured flux point based on the initial individual 1-sigma uncertainty. From 
these spectra we computed the mean values and related errors of the listed parameters. 
Especially at the short wavelength end it shows that not all lines were detected with a 3- 
sigma confidence level. However, since most lines were convincingly measured we decided to 
keep also the tentative detections in our analyses. 

Since we assume LTE the rotational energy levels shall be populated following Boltz- 
mann statistics with a unique temperature for several lines. The involved temperatures T rot 
of the shocked material can be derived from a so-called "rotational diagram". By using the 
results from Eq. [fl plotting the logarithm of N u (J)/(g s gj) against Ej/k (i.e. the formal 
temperature corresponding to the absolute upper energy level of the respective rotational 
transition) and fitting a straight line to the data one finds that the slope of the straight line 
is proportional to — l/T rot . Here, g s denotes the spin degeneracy of each energy level (1 for 
even J, 3 for odd J) and gj = 2 J + 1 is the rotational degeneracy. These numbers assume 
that the ortho-to-para ratio of the involved hydrogen is close to its LTE value of 3 at T rot . A 
deviation from this LTE assumption (i.e., ortho-to-para < 3) would result in a downwards 
displacement of the data points with an odd J-number (ortho-H 2 ) relative to the points with 
even J-numbers (para-H 2 ) and thus creat e a "zig-zag" pattern in the rotational plot (e.g., 



Neufeld et al.l 12009 : IWikenbus et al.ll2000h 



In Fig. [HI [91 and [10] we show rotational plots for the measured lines near RNO IB and 
RNO 1C. In all plots the errors for the data points denote the 1-sigma uncertainty in the 
measured column density. 

The rotational diagram of RNO IB in Fig. [8] shows a clear curvature in the data points. 
This departure from a single straight line is indicative of several temperature components in 
the shocked material. We thus fitted the data with two superimposed temperature regimes 
(hot and warm) that rather represent the minimum and maximum of the involved tem- 
peratures than any distinct intermediate value. The dash-dotted line fits the high energy 
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regime and represents the hot component with a temperature of 2991±596 K. The dotted 
line corresponds to a temperature of 1071±121 K designating the warm component. The un- 
certainties in the temperature denote the 1-sigma confidence level of the fits. As the derived 
temperature is extremely sensitive to the slope of the fitted straight line the corresponding 
errors are rather large. 

In addition to the temperatures one can also derive the total H 2 column density of the 
observed shock within the given aperture. From the y-intersects of the fits we estimate the 
hot component to have iV-f^^H^)^ 1.2 x 10 17 cm -2 and for the warm component we find 
A^ arm (H 2 )^ 5.6 x 10 18 cm" 2 . 

For the column densities derived from the spectrum close to RNO 1C we plotted two ro- 
tational diagrams shown in Fig. [H]and[TDJ First, we fitted the data with a single temperature 
component of 1754±321 K. However, the shocked material is apparently not in LTE as the 
data points with odd and even J-numbers are difficult to fit simultaneously with a straight 
line or a curve. As mentioned above this displacement is indicative of a departure from the 
LTE ortho-to-para ratio of 3. Fig. [TD] shows the same data if one assumes an ortho-to-para 
ratio of 1. Now, a two component fit, similar to that in Fig. [SJ is possible. For the hot com- 
ponent we derive a temperature of 2339±468 K, while the warm component is significantly 
colder with 466±264 K. The total H2 column densities amount to A r | 1 ( 9 t (H2)~ 5.5 x 10 17 cm -2 
and N^ m (R 2 )~ 9.7 x 10 19 cm" 2 , respectively. 

Table H] summarizes the temperatures and total column densities of the different com- 
ponents for both spectra. 



Discussion 



The IRAC images reveal for the first ti me the embedded MIR point source associated 
with IRAS 00338+6312. However, already IWeintraub fc Kastnerl (119931 ) concluded from 
polarimetric observations that an additional source should be present close to the IRAS 
point source position. Furthermore, many authors analyzing the bipolar outflow related to 
this region concluded that a deeply embedded object was m ost likely the driving source 



(lYang et al. 



1995 



Xueta 



objec ts ( jMcMuldroch et al 



2006; IWeintraub & Kastnerl Il993l ) instead of o ne of the RNO 



1995T ) . Previous ground based MIR observations (IPolomski et al. 



20051 ) were possibly not sensitive enough to detect IRAS 00338+6312. 



The observed properties of IRAS 00338+6312 convincingly classify the object as a young 
protostar and sup port the idea tha t it is indeed driving the molecular outflow. From their 
CO measurements Ku et al.l (120061 ) derived a total mass for this bipolar outflow of 1.4 M , 
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indicating that IRAS 00338+6312 is an intermediate-mass protostar. The detection of the 
IRAS source solves also the problem that up to now the SE Ds of RNO IB and RNO 1C 
were unusually steep longward of 25 /im for FUOR objects (IWeintraub fc Kastnerl Il993l ). 
The images presented here have sufficient sensitivity and spatial resolution to show that the 
newly revealed source associated with IRAS 00338+6312 is most probably responsible for 
this flux excess. 



In addition to the IRAS object IWeintraub fc Kastnerl (119931 ) also suspected another 
embedded source which we identify with RNO 1G. Furthermore, the detection of RNO IF 
(IWeintraub et al.lll996l ). which was initiall y assumed to be a density e nhancement rather 
than an independent self-luminous source (IWeintraub fc Kastnerl Il993l ). is confirmed with 
our IRA C images. Only RNO I D, which was thought to be located between RNO IB and 
RNO 1C (jStaude &: Neckellll99lh. does not show u p in our images. This source was, however, 



also not detected by 



Weintraub et all fll996j ). 



Weintraub fc Kastnerl (119931 ) and only tentatively seen in the data from 



It is not clear whether all detected objects are physically confined to the RNO 1B/1C 
region. Especially the membership of the newly discovered objects RNOl IRAC1, RNOl 
IRAC2, and RNOl IRAC3 needs to be confirmed. However, the IRAC observations suggest 
that RNO IB and RNO 1C belong to a small cluster of (partly very) young objects. To 
our knowledg e, although some FUOR s are known or suspected to be in binary or multiple 
systems (e.g.. lReipurth fc Aspinll2004l ). the existence of FUORs in a cluster-like environment 
is a so far un ique finding. Only the F UOR candidate VI 184 Tau (CB34) belongs also to a 



small cluster (IKhanzadvan et al. 



for the RNO objects (jYun et al. 



2002) but its classification as FUOR is far less certain than 



19971 ). In any case the still deeply embedded protostars 



in the close vicinity of RNO 1B/1C put strong constraints on the age of the objects if one 
assumes coeval evolution. 

The magnitudes we derive for the two FUO R objects at 3 . 6 um can be compared to 



previo us ground based measurements at 3.8 /im by lKenyon et al.l (119931 ) and lPolomski et al. 
( 120051 ). However, one certainly has to keep in mind that the IRAC filter is not only different 



in terms of central wavelen gth but also t h e spec tral width differs from that of the ground 
based instruments. While iKenyon et al.l (119931 ) found L^ jm = 6.36 mag and L 3 8flia = 
6.5 mag for RNO IB and RNO 1C, respectively, IPolomski et all (bood ) observed RNO 1C 
to be slightly brighter than RNO IB (L 3 , 8 ^ m = 6.27 mag vs. L 3 , 8 ^m = 6.42 mag). Also in 
our measurements RNO 1C seems to be the bright er component. Howe ver, both objects 
appear to be fainter in our observations compared to lPolomski et al.l (120051 ). Apart from the 
differences in the filter properties intrinsic variations in the luminosity of the FUORs and/or 
changing local extinction effects can account for the apparent variability in these objects. 



- 11 - 



The detection of H2 rotational line e mission can in general be attributed to collisiona l 



excitation from C- or J-type shocks (e.g.. iDraine &: McKedll993l ; iMoro-Martm et al.ll200ll ). 



C-type shocks are magnetic and less violent towards molecules in comparison to the mostly 
hydrodynamic J-type shocks that can dissociate H 2 molecules even at lower shock velocities. 
From our rotational diagrams we derived shock temp eratures that are similar to those found 
in the outflow from Cepheus A (jFroebrich et al.ll2002l ). These authors could fit a two compo- 
nent C-shock model to their data and derived shock velocities between 25 and 30 kms -1 for 
a cold and a hot component, respectively. By comp aring our results dire ctly to theoretical 
shock models we find that the C-shock Model 1 from iTimmermannl (119981 ) can explain both 
the hot and warm component in the spectrum close to RNO IB. Considering our derived 
temperatures as an upper and lower limit for the sh ocked material , shock velocities between 
15 and 30 kms -1 are required. However, Model 1 of ITimmermannl (119981 ) predicts on ortho- 
to-para ratio of ~2 at the low velocity/low te mperature limit which is not directly evident 
from our observations. More recent models by IWilgenbus et al.l (120001 ) predi ct lower shock 
tempe ratures in C-type shocks for the same shock velocities in comparison to iTimmermann 
(119981 ) . Following their models, our hot component in FigJHJ requires velocities exceeding 
40 kms -1 . Slightly slower velocities (~35 kms" 1 ) are already needed to explain the hot 
temperature component in the shock close to RNO 1C (Fig. fTOl) . The warm component 
in this diagram represents shock velocities between 15 kms -1 and 20 kms -1 . Although it 
is clear that some uncertainties between the observations and theory still remain, C-type 
shocks seem to provide a solid explanation for the observed shock temperatures. 

The measurement of an ortho-to-para ratio smaller than 3 in Fig. [TD] implies that despite 
the high temperatures the gas has not yet reached equilibrium between the ortho and para 
states. Thus, the observed ratio is the legacy of the temperature history of the gas (see, e.g., 
Neufeld et al.l Il998l . 120061 ) and transient heating by a recently passing shock wave already 
caused higher temperatures. As the corresponding region on the sky lies closer to the IRAS 
source than the region probed in Fig. [HI the observations suggest that we see the signatures 
of at least two shock waves: The first shock wave is probed close to RNO IB (Fig. [HJ where 
we observe several temperature components that have apparently reached their equilibrium 
state. A more recent shock wave is seen close to RNO 1C (Fig. ITUT) where the ratio of the 
line column densities hints also towards high shock velocities and temperatures but the gas 
is not yet in thermal equilibrium. 



5. Conclusions and Future Prospects 

Our conclusions can be summarized as follows: 
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• We detected and resolve for the first time the MIR point source associated with IRAS 
00338+6312 which appears to be an embedded intermediate-mass protostar driving 
the known molecular outflow in the RNO 1B/1C region. 

• The detection of additional (partly previously unknown) point sources suggests that the 
FUOR objects RNO 1B/1C belong to a young small stellar cluster. To our knowledge, 
RNO 1B/1C are the only well-studied and confirmed FUORs that apparently belong 
to a cluster-like environment. 

• All but two objects were detected in all four IRAC bands and their colors are consistent 
with Class O/I-II objects. The two objects that were not detected at 3.6 /im (including 
IRAS 00338+6312) are still very deeply embedded protostars. 

• Having apparently extremely young objects in the direct vicinity of the FUORs con- 
firms the suspected young age for this type of objects although their MIR colors are 
consistent with Class II objects. 

• The two MIR spectra of the region bear clear signs of a dense icy and dusty circum- 
stellar environment as solid state features are seen in absorption. 

• The spectra show also H 2 emission lines from purely rotational transitions. We presume 
that these lines arise from shocked material within the molecular outflow. The derived 
shock temperatures and velocities are in agreement with C-type shock models. 

• The observations of the H 2 lines suggest that the outflow lies in front of RNO 1B/1C 
as otherwise the high optical depth towards these objects would have prevented the 
detection. 

• While in one spectrum the gas probed by the H 2 line emission seems to be in LTE, 
the other spectrum shows a deviation from the expected LTE ortho-to-para H 2 ratio. 
This indicates the presence of at least two shock waves, the most recent one responsible 
for the non-LTE line ratios. 



The results presented here advantageously combine space-based MIR imaging and spec- 
troscopy. While the images give insights into the photometric properties of several young 
embedded objects and also FU Orionis type stars the spectra provide information on the 
dusty and icy circumstellar environment of the young cluster as well as on shocked gas 
within a molecular outflow. In this context the existence of FUORs in a cluster-like en- 
vironment needs to be pointed out. It is still debated whether all young low-mass stars 
undergo an FUOR phase o r whether these objects form a special sub-group of YSOs (e.g., 



Hartmann &: Kenyonlll996l ). Our results suggests that the FUOR phenomenon occurs also 
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in young (small) stellar clusters making it thus possibly more common than so far expected. 
In addition, the co-existence of FUORs and very young protostars in the same environment 
strengthens the idea that the FUOR-phase of YSOs is linked to the early stages of the star 
formation process. 

Possible future investigations of the RNO 1B/1C region should address whether the 
observed objects do indeed belong to the same region and whether all of them are young stars. 
Furthermore, a detailed study of the shocked material including high-resolution narrow band 
imaging and spectroscopic mapping (e.g., with Spitzer IRS) will provide deeper insights into 
its extension, physical conditions, and its relation to the molecular outflow. A complete 
analysis of the ice and dust features in the observed spectra is currently underway and will 
be presented in an upcoming paper (Quanz et al. in preparation). 

S. P. Quanz kindly acknowledges support from the German Friedrich-Ebert-Stiftung. 
We thank H. Beuther, R. Mundt, A. Carmona and S. Birkmann for useful discussions and 
conversations. An anonymous referee provided valuable criticism that helped improve the 
paper. This work is based on observations made with the Spitzer Space Telescope, which 
is operated by the Jet Propulsion Laboratory, California Institute of Technology under a 
contract with NASA. This research has made use of the SIMBAD database, operated at 
CDS, Strasbourg, France. 

Facilities: Spitzer 



A. The non-detection of the 28.22 /im emission line 

In the following we analyze whether the non-detection of the 28.22 /im in either spectrum 
is in agreement with the models. The two component models in Fig. [H] and (TU] can be used 
to predict the column densities of the S(0) emission line. One expects 

, N S(0) r - 2 , 10 „ _ 2 

log — — cm « 18.3 cm 

gjgs 

log— ^[cnr 2 l « 18.8 cm" 2 

gjgs 

for the spectrum close to RNO IB and RNO 1C, respectively. With g s =l and gj—5 and 
applying Eq. (JTJ the following line intensities are derived: 

I(2)ib « 1.64 • 10~ 6 erg s" 1 cm' 2 sr" 1 

I(2)ic « 5.20 ■ 10" 6 erg s" 1 cm" 2 sr" 1 



-14- 



The ^4-coefficient for the transition is 2.94 ■ 10 s (IWolniewicz et al.lll998l ). 



Taking into account the aperture size for the long wavelength high-resolution module 
of the spectrograph of 247.53 arcsec 2 (~ 5.82 ■ 10~ 9 sr) the predicted integrated line fluxes 
amount to: 

Fib ~ 9.57- 10~ 22 Wcm~ 2 
r21 Wcm- 2 



Fic « 3.02 • 10" 



To estimate the peak of the 28.22 /im line in Jansky we assume that the spectral resolu- 
tion in the high-resolution modules is constant and that the FWHM of the S(0) line can be 
extrapolated from the FWHM of the S(l) line at 17.03 pm. Since we fit a Gaussian profile 
to the emission lines and we want to derive the peak flux of this Gaussian, we also have to 
take into account the relation between the FWHM we measure for the line and the a of the 
profile. This relation is FWHW=v / 8 ln2 ■ a. With the expected FWHM of 

Ai/ib « 2.84- 10 10 Hz 

Az/ 1G « 2.12 • 10 10 Hz 
we find thus corresponding peak fluxes of: 

F% ak « 0.034 Jy 

F^ ak ^0.152 Jy 



These values have to be compared to the measured uncertainties in the spectra. The 
mean 1—a level in the spectral range between 28.0 and 28.4 /zm is, however, 

a 1B = 0.055 

ctic = 0.082 

so that even in the more favorable case of RNO 1C the emission line is not expected to be 
detected with a confidence level greater than ^1.8 cr. 
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Table 1. Journal of observations. 



Instrument 


RA (2000) a 


DEC (2000) a 


AOR Key 


Filter or Module 


Ramp duration / 

-44- *~iT 1 

jf 1 ui oycies 


Frame Time / 

-44- /~\T HVQnioc / 

jf 1 oi n anies / 
Dither Positions 


Date 


Spitzer IRAC 


00:36:45.8 


+63:28:56 


5027072 


3.6, 4.5, 5.8, 8.0 //m 




0.4/1/4 


2003-12-20 


Spitzer IRS 


00:36:46.34 b 


+63:28:53.76 b 


6586624 


Low Resolution 


6 sec / 3 




2004-01-07 




00:36:46.89 c 


+63:28:58.44 c 












Spitzer IRS 


00:36:46.34 b 


+63:28:53.76 b 


6586624 


High Resolution 


6 sec / 5 




2004-01-07 




00:36:46.89 c 


+63:28:58.44 c 













a Average slit position of low-resolution spectrograph. 
b close to RNO IB 
c close to RNO 1C 



Table 2. Apparent brightness at different wavelengths for objects that were detected in at least three IRAC bands. 
The coordinates are measured in the 3.6 /xm image. Only the positions of the IRAS source and RNO 1G were 

measured in the 4.5 /xm exposure. 



No. 


Object name 


RA (J2000) 


DEC (J2000) 


3.6 /iin 


4.5 /im 


5.8 yum 


8.0 fim 










[mag] 


[mag] 


[mag] 


[mag] 


1 


RNO IB 


00:36:46.05 


+63:28:53.29 


7.16+0.13 


6.67+0.06 


5.76+0.12 


5.01+0.09 


2 


RNO 1C 


00:36:46.65 


+63:28:57.90 


6.56+0.01 


6.04+0.01 


5.58+0.01 


4.61+0.02 


3 


RNO IF 


00:36:45.74 


+63:29:04.09 


10.28+0.09 


9.46+0.08 


8.58+0.05 


8.20+0.03 


4 


RNO lG a 


00:36:47.14 


+63:28:49.95 




10.33+0.06 


8.74+0.04 


8.05+0.07 


5 


IRAS 00338+6312 


00:36:47.34 


+63:29:01.61 




9.05+0.07 


7.19+0.05 


6.72+0.03 


6 


RNOl IRAC1 


00:36:48.44 


+63:28:39.98 


13.72+0.15 


11.82+0.14 


10.93+0.20 


9.65+0.04 


7 


RNOl IRAC2 


00:36:47.90 


+63:28:36.30 


12.00+0.08 


10.86+0.10 


10.86+0.16 


9.82+0.14 


8 


RNOl IRAC3 


00:36:47.85 


+63:28:41.23 


13.07+0.10 


11.78+0.08 


10.74+ 0.10 


9.74+0.05 



Embedded YSO from IWeintraub fc Kastnerl (119931 ) 



Table 3. Properties of the observed Hydrogen emission lines close to RNO IB (upper part) and RNO 1C (lower 
part). The numbers in parentheses denote powers of ten. The 12.28 jum emission line close to RNO IB was measured 
in the low- and high-resolution part of the spectrograph. See text for discussion about the apparent flux difference 

between the two measurements. 



Transition 


Wavelength 


bnergy 


A-L/oemcient 


Beam Size 


Line Flux d 


T * T -i- ' X- d 

Line Intensity 


Oolumn Density ' 








Ej/k [K] 


[s- 1 ] 


[arcsec 2 ] 


[W cm' 2 ] 


[erg cm -2 s _1 sr _1 ] 


[cm -2 ] 










Ho 


lines close to RNO IB 






S(l) J= 


=3-1 


1 7 0348 


1 01 Fi 08 




53.11 


7.91(-21)±3.8(-22) 


a 047 c\_i_o o(_f.\ 


1 44(1 Ql+fi 8(1 71 

X .tit: I ± U l _I_VJ .Oil 1 J 


S(2) J= 


=4-2 


12 2786 


1681 63 


2.75(-9) 


53.11 


1.33(-20)±5.2(-22) 


1 06 (-41 ±4 l(-6l 


3 00(181+1 2(171 












39.96 


6.75(-21)±1.2(-21) 


7 18(-5l+1 3(-5l 

1 . 1UI KJ J _1 — L.ul U j 


o 03d 81+3 7(171 


S(3) J= 


=5-3 


9.6649 


2503 73 


9 83(-9l 


39.96 


2.00(-20)±2.7(-21) 


2 13 Ml ±2 9 (-5) 


1 33(181+1 8(171 


S(4) J= 


=6-4 


8.0251 


3474.48 


2.64(-8) 


39.96 


9.12(-21)±1.9(-21) 


9.71(-5)±2.1(-5) 


1.87(17)±4.0(16) 


S(5) J= 


=7-5 


6.9095 


4585.94 


5.88(-8) 


39.96 


2.13(-20)±2.6(-21) 


2.26(-4)±2.8(-5) 


1.68(17)±2.1(16) 


S(6) J= 


=8-6 


6.1086 


5829.66 


1.14(-7) 


39.96 


7.04(-21)±1.3(-21) 


7.50(-5)±1.4(-5) 


2.54(16)±4.8(15) 


S(7) J= 


=9-7 


5.5112 


7196.20 


2.00(-7) 


39.96 


1.32(-20)±6.1(-21) 


1.41(-4)±6.5(-5) 


2.46(16)±1.1(16) 










H 2 


lines close to RNO 1C 






S(l) J= 


=3-1 


17.0348 


1015.08 


4.76(-10) 


53.11 


3.26(-21)±2.7(-22) 


2.61(-5)±2.1(-6) 


5.92(18)±4.8(17) 


S(2) J= 


=4-2 


12.2786 


1681.63 


2.75(-9) 


53.11 


5.17(-21)±7.7(-22) 


4.14(-5)±6.2(-6) 


1.17(18)±1.7(17) 


S(3) J= 


=5-3 


9.6649 


2503.73 


9.83(-9) 


39.96 


7.52(-21)±2.0(-21) 


8.00(-5)±2.1(-5) 


4.98(17)±1.3(17) 


S(4) J= 


=6-4 


8.0251 


3474.48 


2.64(-8) 


39.96 


1.34(-20)±4.1(-21) 


1.43(-4)±4.4(-5) 


2.74(17)±8.5(16) 


S(5) J= 


=7-5 


6.9095 


4585.94 


5.88(-8) 


39.96 


1.10(-20)±3.3(-21) 


1.17(-4)±3.5(-5) 


8.69(16)±2.6(16) 



1 Energy level of the upper state following iJennings et al.l (119871 ) 



Wolniewicz et al.l (119981 ) 

Corresponds to the instrument aperture size in case of the short high-resolution module (53.11 arcsec 2 , see Spitzer Handbook). 
For all other transitions (i.e., short low-resolution module) the number denotes the extracted beam size during data reduction process 
(39.96 arcsec 2 ). 

d Corrected for extinction using Ay = 3.55 mag. 



c Upper energy state 



Table 4. Upper and lower limits on the temperatures, and the total H 2 column densities of the observed shocked 
components. For the data derived close to RNO 1C the results from the two component fit and the single component 

fit with different ortho-to-para ratios are given. 



Component 


T [K] lo 


gio N(H 2 ) [cm" 2 ] 


RNO IB 


Hot a 


2991±596 


17.09±0.36 


Warm a 


1071±121 


18.75±0.09 




RNO 1C 




Hot b 


2339±468 


17.74±0.31 


Warm b 


466±264 


19.99±1.17 


Single a 


1754±321 


17.98±0.31 



a ortho-to-para ratio of 3 
b ortho-to-para ratio of 1 
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Fig. 1. — The slits of the different modules of the Spitzer/IRS overplotted on the inverted 
2MASS Ks-filter image. The two FUOR objects RNO IB and RNO 1C are not centered in 
the slits (north is up, east to the left). 



-23- 




Fig. 2.— The RNO 1B/1C region in the 2MASS Ks-filter (left) and the 5.8 /an IRAC filter 
(right). The images are inverted and the contrast is set to make also fainter objects visible. 
The IRAC image contains more sources and especially the IRAS source clearly shows up at 
5.8 /zm while it seems hidden behind a dust lane in the 2MASS image. The magnitudes of 
the individual sources are listed in Table [2j 
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IRAC color— color plot 
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RAC color-color plot 



5 - 



T 3 

2 



.-2-4-.r 
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0.0 



0.5 1.0 1.5 

[4.5]-[5.8] 



2.0 



Fig. 3. — IRAC color-color plots for the objects listed in Tabl e [2 



the regions corresponding to Class II objects in lHartmann et al. 



The d ashed boxes indicate 
J2005h. The dashed - dotte d 

region in the right plot defines the positio n of Class 0/1 sources in lHartmann et al. ( 2005 ). 
In the left plot the Class 0/1 objects from lHartmann et al.l (120051 ) have the same [5.8]-[8.0] 
colors as the Class II objects but redder [3.6]- [4.5] colors. Thus, they lie "above" the dashed 
box. The reddening vectors correspond to Ay = 30 mag and are based on a Vega-like 
spectrum and the reddening law given bv iMathisI Jl99ph . 
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RNQ f'B with SpHm/lRS 




7 E 9 10 11 >2 13 
RNO 1C with 5pi:zcr/'RS 



12 13 14 15 16 17 13 19 20 



20 22 



24 26 2£ 
Wovelengtli 



30 32 34 




5 6 7 6 9 10 11 ■>■! 13 
Wavderqth [^"i] 



12 13 14 15 16 17 18 19 20 
Wavelength [)irn] 




20 22 



Fig. 4. — Spitzer/IRS spectra close to RNO IB (upper row) and RNO 1C (lower row). 
The leftmost column shows the low-resolution data while the other two columns present the 
high-resolution spectra. The positions of the H 2 -emission lines are indicated by the arrows. 
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Fig. 5. — Comparison of the optical depth of the 10 /im absorption feature with a typical 
ISM absorption feature. The solid lines show the optical depth observed towards the RNO 
comple x while the dashed lines represent the (scaled) optical depth towards the galactic 
center (IKemper et al.l I2004J ) . The absorption in the RNO region can not be solely caused 
by pure ISM silicate grains as both absorption features seem to be slightly shifted towards 
shorter wavelengths. Furthermore, while the spectrum close to RNO 1C shows additional 
absorption longwards of 11 /im (possibly blended with H 2 ice at 13/iin), the feature close 
to RNO IB is too narrow and an additional emission feature might contribute longwards of 
~10.5 /im. 
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-0.02 f , . , Lj 

6.639 6.911 

Wovelength [^m] 



0.03 
0.02 
0.01 
0.00 



-0.01 
1 1.739 



S(2) Line 



— ^/^ 




1 





12.298 
Wovelength [/irr 



S(6) Line 



-0.005 r , ,__L 

5.814 6.120 

Wovelength [jj.fr 




0.05 
0.04 
0.03 
0.02 
0.01 
0.00 
-0.01 
9.102 



0.025 
0.020 



S(3) Line 



S(4) Line 







i 



9.670 
Wovelength [^m] 



S(7) Line 




-0.005 Lj ill ^ i J 

5.282 5.529 5.776 

Wovelength [/^ni] 




8.023 
Wavelength [,iim] 



Fig. 6. — Gaussian fits to the H2 lines in the spectrum close to RNO IB after continuum 
subtraction. The line at 12.28 fim shown here was measured in the low-resolution spectrum. 
Table[3j however, contains the flux of this line measured in both, the low- and high-resolution 
part of the spectrum. The line at 6.11 /im is blended with the ice feature in this region. 
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S(1) Line S(2) Line S(3) Line S(4) Line 




16.875 17.039 17.203 12.242 12.277 12.312 9.289 9.652 10.015 7.435 8.002 8.570 

Wavelength [/im] Wavelength [/xm] Wavelength [^m] Wavelength [/im] 



S(5) Line 
0.06 1 1 1 I 1 




-0.02 1 

6.747 6.896 7.045 

Wavelength [jLim] 



Fig. 7. — Same as Fig. [6] but for the spectrum close to RNO 1C. Here, the 12.28 /zm line was 
only detected in the high-resolution spectrum. 
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Fig. 8. — Rotational diagram of the H2 lines shown in Fig. El The data are well constrained 
by the superposition of a hot and a warm component providing an upper and a lower limit 
for the temperature range of the shocked gas. The diamond at ^1682 K corresponds to the 
12.28 /im line measured in the high-resolution module. The observed line flux, and hence the 
derived column density, is slightly higher than that in the low-resolution spectrum. For the 
fit, however, the low-resolution data point was used as all other higher energy data points 
were also obtained with the low-resolution spectrograph and uncertainties related to the 
different apertures and their orientations for the high- and low-resolution modules are thus 
minimized. 
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Fig. 9. — Rotational diagram of the H2 lines shown in Fig. El The apparent zig-zag pattern 
is indicative of an ortho-to-para ratio smaller than 3. 
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Fig. 10. — Same as Fig.[9]but now for an assumed ortho-to-para ratio of 1. A two component 
fit similar to that in Fig. [H] is now possible. 
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RNO 1 B/C complex seen with Spitzer/IRAC 




RIGHT ASCENSION 
CENTRE: R.A. 00 i, 36 m 4e*46 DEC +G3°28'54: , 9 (J20O0) 




Fig. 11. — Color composite of the RNO 1B/1C region based on the three IRAC bands at 
3.6, 4.5, and 8.0 /im. The image is centered on RNO 1C. We applied a logarithmic color 
stretch to enhance the fainter sources. The extended emission east of the IRAS source which 
is predominantly seen in the green 4.5 /zm filter is possibly related to H 2 emission within the 
bipolar outflow of this object. 



